Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

Profili di espressione genetica ed epigenetica ad impatto prognostico e predittivo nel carcinoma squamoso e nelle lesioni potenzialmente maligne del cavo orale

Il carcinoma squamoso orale (CSO) �� spesso preceduto da lesioni definite potenzialmente maligne tra cui la leucoplachia e il lichen ma una diagnosi precoce avviene ancora oggi in meno della met�� dei casi. Inoltre spesso un paziente trattato per CSO svilupper�� secondi tumori. Scopo del lavoro di ricerca �� stato: 1) Studiare, mediante metodica di next generation sequencing, lo stato di metilazione di un gruppo di geni a partire da prelievi brushing del cavo orale al fine di identificare CSO o lesioni ad alto rischio di trasformazione maligna. 2) Valurare la relazione esistente tra sovraespressione di p16INK4A e presenza di HPV in 35 pazienti affetti da lichen 3) Valutare la presenza di marker istopatologici predittivi di comparsa di seconde manifestazioni tumorali 4) valutare la relazione clonale tra tumore primitivo e metastasi linfonodale in 8 pazienti mediante 2 metodiche di clonalit�� differenti: l���analisi di mtDNA e delle mutazioni del gene TP53. I risultati hanno mostrato: 1) i geni ZAP70 e GP1BB hanno presentato un alterato stato di metilazione rispettivamente nel 100% e nel 90,9% di CSO e lesioni ad alto rischio, mentre non sono risultati metilati nei controlli sani; ipotizzando un ruolo come potenziali marcatori per la diagnosi precoce nel CSO. 2)Una sovraespressione di p16INK4A �� risultata in 26/35 pazienti affetti da lichen ma HPV-DNA �� stato identificato in soli 4 campioni. Nessuna relazione sembra essere tra sovraespressione di p16INK4A e virus HPV. 3)L���invasione perineurale �� risultato un marker predittivo della comparsa di recidiva locale e metastasi linfonodale, mentre lo stato dei margini chirurgici si �� rilevato un fattore predittivo per la comparsa di secondi tumori primitivi 4) Un totale accordo nei risultati c����� stato tra analisi di mtDNA e analisi di TP53 e le due metodiche hanno identificato la presenza di 4 metastasi linfonodali non clonalmente correlate al tumore primitivo.

Evidenza del coinvolgimento di geni nell'insorgenza delle schisi labio-palatine

The CL/P are the most common and easily recognizable craniofacial malformations with a complex etiology that requires the involvement of genetic and environmental components. The analysis of the genetic component shows more than 14 loci and genes involved in the onset of the disease. I���ve selected and investigated some of the possible candidate genes for CL/P. MYH14 gene, that maps on chromosome 19, on the OFC3 locus, and shows a strong homology with MYH9 gene. I���ve also investigated TP63 and MID1 genes, that are responsible respectively for EEC syndrome and Opitz syndrome, both of them presenting cleft. I���ve also decided to investigate JAG2 because TP63 product regulates the this gene, and both of them are component of the Notch signalling pathway. I���ve, also, studied the MKX and LMO4 genes. MKX is an important development regulator that is highly expressed in palatal mesenchyme, and map in the region responsible for Twirler mutation that cause cleft in mouse. LMO4 is necessary for neural tube development and cooperating with Grhl3, promotes cellular migration during morphogenetic events like ���in utero��� cleft healing. Low folate levels and high levels of homocysteine increase the risk of cleft, genes involved in their metabolism may be of interest in cleft occurrence. I���ve decided to investigate BHMT and CBS genes coding for enzymes involved in homocysteine metabolism. I���ve also investigated BHMT2 gene that maps close to BHMT and presents with him a 73% of homology. I���ve performed a linkage analysis using SNPs mapping in the genes and their boundaries, for each gene, for MKX and LMO4 I���ve also performed a sequencing analysis. My results for MID1 and CBS genes support the hypothesis of a possible role of these genes in cleft. I���ve found borderline association values for JAG2, MKX and LMO4 genes.

Immunohistochemical and Molecular Prognostic/Predictive Markers in Neoplastic Diseases

Traditional morphological examinations are not anymore sufficient for a complete evaluation of tumoral tissue and the use of neoplastic markers is of utmost importance. Neoplastic markers can be classified in: diagnostic, prognostic and predictive markers. Three markers were analyzed. 1) Insulin-like growth factor binding protein 2 (IGFBP2) was immunohistochemically examined in prostatic tissues: 40 radical prostatectomies from hormonally untreated patients with their preoperative biopsies, 10 radical prostatectomies from patients under complete androgen ablation before surgery and 10 simple prostatectomies from patients with bladder outlet obstruction. Results were compared with ��-methylacyl-CoA racemase (AMACR). IGFBP2 was expressed in the cytoplasm of untreated adenocarcinomas and, to a lesser extent, in HG-PIN; the expression was markedly lower in patients after complete androgen ablation. AMACR was similarly expressed in both adenocarcinoma and HG-PIN, the level being similar in both lesions; the expression was slightly lower in patients after complete androgen ablation. IGFBP2 may be used a diagnostic marker of prostatic adenocarcinomas. 2) Heparan surface proteoglycan immunohistochemical expression was examined in 150 oral squamous cell carcinomas. Follow up information was available in 93 patients (range: 6-34 months, mean: 19��7). After surgery, chemotherapy was performed in 8 patients and radiotherapy in 61 patients. Multivariate and univariate overall survival analyses showed that high expression of syndecan-1 (SYN-1) was associated with a poor prognosis. In patients treated with radiotherapy, such association was higher. SYN-1 is a prognostic marker in oral squamous cell carcinomas; it may also represent a predictive factor for responsiveness to radiotherapy. 3) EGFR was studied in 33 pulmonary adenocarcinomas with traditional DNA sequencing methods and with two mutation-specific antibodies. Overall, the two antibodies had 61.1% sensitivity and 100% specificity in detecting EGFR mutations. EGFR mutation-specific antibodies may represent a predictive marker to identify patients candidate to tyrosine kinase inhibitors therapy.

Role of new molecular approaches for the early diagnosis of bladder cancer in clinical practice

There is an urgent need to improve the performance of urine cytology for the diagnosis of bladder cancer. In preliminary studies, telomerase activity evaluated by telomeric repeat amplification protocol (TRAP) assay and chromosomal aneuploidy detected by fluorescence in situ hybridization (FISH) in the diagnosis of bladder cancer have produced important results. Urine cell-free (UCF) DNA has also been proposed as a potential marker for early bladder cancer diagnosis. In the first study the diagnostic performance of TRAP assay and FISH analysis was assessed, while the second study evaluated the potential role of UCF DNA integrity in early bladder cancer diagnosis. In the first cross-sectional study, 289 consecutive patients who presented with urinary symptoms underwent cystoscopy and cytology evaluation. In the second study, UCF DNA was isolated from 51 bladder cancer patients, 46 symptomatic patients, and 32 healthy volunteers. c-Myc, BCAS1 and HER2 gene sequences longer than 250 bp were quantified by real time PCR to verify UCF DNA integrity. In the first study, sensitivity and specificity were 0.39 and 0.83, respectively, for cytology; 0.66 and 0.72 for TRAP; 0.78 and 0.60 for the cytology and TRAP combination; 0.78 and 0.78 for the cytology, TRAP and FISH combination; and 0.65 and 0.93 for the TRAP and FISH combination. In the second study, at the best cutoff of 0.1 ng/��l, UCF DNA integrity analysis showed a sensitivity of 0.73 and a specificity of 0.84 in healthy individuals and 0.83 in symptomatic patients. The preliminary results suggest that these biomarkers could potentially be used for the early diagnosis of bladder cancer, especially in high-risk populations (e.g, symptomatic individuals exposed to occupational risk) who may benefit from the use of noninvasive diagnostic tests in terms of cost-benefit.

Microbial ecology of biotechnological processes

The investigation of phylogenetic diversity and functionality of complex microbial communities in relation to changes in the environmental conditions represents a major challenge of microbial ecology research. Nowadays, particular attention is paid to microbial communities occurring at environmental sites contaminated by recalcitrant and toxic organic compounds. Extended research has evidenced that such communities evolve some metabolic abilities leading to the partial degradation or complete mineralization of the contaminants. Determination of such biodegradation potential can be the starting point for the development of cost effective biotechnological processes for the bioremediation of contaminated matrices. This work showed how metagenomics-based microbial ecology investigations supported the choice or the development of three different bioremediation strategies. First, PCR-DGGE and PCR-cloning approaches served the molecular characterization of microbial communities enriched through sequential development stages of an aerobic cometabolic process for the treatment of groundwater contaminated by chlorinated aliphatic hydrocarbons inside an immobilized-biomass packed bed bioreactor (PBR). In this case the analyses revealed homogeneous growth and structure of immobilized communities throughout the PBR and the occurrence of dominant microbial phylotypes of the genera Rhodococcus, Comamonas and Acidovorax, which probably drive the biodegradation process. The same molecular approaches were employed to characterize sludge microbial communities selected and enriched during the treatment of municipal wastewater coupled with the production of polyhydroxyalkanoates (PHA). Known PHA-accumulating microorganisms identified were affiliated with the genera Zooglea, Acidovorax and Hydrogenophaga. Finally, the molecular investigation concerned communities of polycyclic aromatic hydrocarbon (PAH) contaminated soil subjected to rhizoremediation with willow roots or fertilization-based treatments. The metabolic ability to biodegrade naphthalene, as a representative model for PAH, was assessed by means of stable isotope probing in combination with high-throughput sequencing analysis. The phylogenetic diversity of microbial populations able to derive carbon from naphthalene was evaluated as a function of the type of treatment.